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Abstract 

We compute atomistically the heat conductance 
for ultra-thin pristine silicon nanowires (SiNWs) 
with diameters ranging from 1 to 5 nm. The room 
temperature thermal conductance is found to be 
highly anisotropic: wires oriented along the (110) 
direction have 50-75% larger conductance than 
wires oriented along the (100) and (111) directions. 
We show that the anisotropics can be qualitatively 
understood and reproduced from the bulk phonon 
band structure. Ab initio density functional the- 
ory (DFT) is used to study the thinnest wires, but 
becomes computationally prohibitive for larger di- 
ameters, where we instead use the Tersoff empiri- 
cal potential model (TEP). For the smallest wires, 
the thermal conductances obtained from DFT- and 
TEP calculations agree within 10%. The pre- 
sented results could be relevant for future phonon- 
engineering of nanowire devices. 



Thermal transport properties of silicon 
nanowires (SiNWs) are crucial for many of 
their future applications. A high thermal con- 
ductance of SiNWs is desirable in nano-electronic 
components, while, on the other hand, a low ther- 
mal conductance is beneficial for thermoelectric 
applications. Indeed, recent experiments have 
demonstrated that SiNWs can have a significantly 
higher thermoelectric performance than bulk sili- 
con [TJ [5] , due to a reduced thermal conductance 



of the nanowires [3]. Raman microscopy studies 
have revealed significant phonon confinement 
effects in SiNWs down to 4 nm in diameter [3]. 
Very thin SiNWs with diameters below 5 nm have 
been synthesized by several groups [31 El El E] 
typically oriented in the (100), (110), (111), or 
(112) directions. 

Most theoretical works on phonon transport in 
nanowires have been based on continuous elastic- 
ity theory and the Boltzmann transport equation, 
e.g. [SI [5]. Recently, however, some theoretical 
work based on atomistic models has been reported: 
Mingo et al. obtained good agreement with exper- 
iments for wires with diameters > 35 nm [TU1 \H\ ; 
Wang et al. [12] studied the diameter dependence 
of pristine (100) wires, and Zhang et al. [13] an- 
alyzed the interface scattering between nanowires 
and bulk contacts. The effects of an amorphous 
coating on the transmission was studied in Ref. 

An important design parameter for SiNW de- 
vices is the crystal orientation of the wires, and 
it is thus highly relevant to compare the thermal 
transport in the different wire directions. Sev- 
eral recent works have already addressed the di- 
rectional dependence on the electronic transport 
properties [13 [H El [U QJ]. To our knowl- 
edge, no experimental comparison of heat conduc- 
tion in different wire orientations in SiNWs has 
been published. In wires connected to bulk con- 
tacts, Zhang et al. [13] found only insignificant dif- 
ference between (100) and (111) wires. Early mea- 



1 



surements [20] [2T] of anisotropic heat conduction 
and phonon focusing in bulk silicon was analyzed 
theoretically [22J for boundary-scattered phonons 
in macroscopic samples. However, in very narrow 
SiNWs with smooth surfaces there will be no sur- 
face scattering, since the surface structure is incor- 
porated in the phonon modes. In the ideal ballistic 
limit, all phonon modes will propagate through the 
wire with a transmission probability of one. 

In this work we compute the thermal conduc- 
tance of ultra thin silicon nanowires (SiNWs) with 
diameters from 1-5 nm oriented along the (100), 
(110), and (111) directions. As a first step we 
use both DFT calculations [53] and empirical po- 
tentials [5H [55J to calculate the dynamical ma- 
trix. Next, the Landauer transmission and ther- 
mal conductance are calculated. We find that 
the room temperature thermal transport is highly 
anisotropic with (110) wires having up to two times 
larger thermal conductance than (100) and (111). 
Finally, we show that the anisotropics can be un- 
derstood based on the underlying bulk phonon 
band structure. 

We use atomic orbital DFT calculations [53J to 
study the smallest wires (diameter ~ 1 nm) in the 
(100) and (110) directions. After an initial relax- 
ation, each atom, /, is displaced by Qj^ in direc- 
tion [i = {x,y, z} to obtain the forces, Fj u (Qj fi ) 1 
on atom J ^= I in direction v |31j . The dynamical 
matrix, K, is then found by finite differences [32] 

K = d 2 E = F Ju (Q I ^-Fj,(-Q I ^ 

(1) 

with E being the total energy. The intra-atomic 
elements are calculated by imposing momentum 
conservation, such that Kj^j^ = — J2k^i Kin,Kv 
The forces are relatively long ranged and we found 
it necessary to use periodically repeated super cells 
containing at least 5 unit cells, corresponding to 
a minimum length L$c > 19 A in the (110) di- 
rection. This implies that the calculations become 
very time consuming and puts strong limitations 
on the diameter range to be investigated. 

Since the large-diameter DFT calculations be- 
come computationally prohibitive, one needs a 
faster computational scheme. To this end we 
use the Tersoff empirical potential model (TEP) 
[54"l 133] as implemented in the "General Lattice 
Utility Program" (gulp) [25]. We use GULP to 



relax the atomic structure and to directly out- 
put the dynamical matrix, K, for the relaxed sys- 
tem. The interactions in the TEP are highly local, 
and the calculations can be carried out on periodi- 
cally repeated super cells containing only two unit 
cells. The TEP calculations are much faster (> 10 4 
times) and much less memory demanding than the 
DFT calculations. In order to validate the empiri- 
cal scheme, we present below comparisons between 
the DFT and TEP schemes. 

Experimentally fabricated SiNWs are typically 
passivated with either amorphous silicon oxide, 
SiC>2, [7J or with hydrogen 5J. The influence of 
amorphous SiC>2 passivation on phonon transport 
was investigated by Mingo and Yang 14J. However, 
atomistic models of amorphous SiC>2 are demand- 
ing and often one uses hydrogen in calculations of 
the electronic transport to passivate the silicon dan- 
gling bonds [26] [27] [28]. The surface passivation 
is crucial for the electronic properties since the Si 
dangling bonds form states in the bandgap possibly 
leading to metallic SiNWs [29 . Passivation by hy- 
drogen suppresses surface reconstruction and leads 
to semiconducting wires. We thus include hydro- 
gen when we model SiNWs using DFT calculations. 
We show below, that it is not necessary to include 
H-passivation in the TEP model. 

We formally divide our system into a left, central, 
and right region. The left and right contacts are 
modeled as two semi-infinite pristine wires. In the 
present study, the central region is also a pristine 
wire, but it could in principle include any struc- 
tural defects [H] [3H [35] as we shall describe in 
a forthcoming paper [36j . In the limit of a small 
temperature difference between the left and right 
regions, the thermal conductance is 34, |37J [35] 

fr2 i-oo e huj/k B T 

G(T) = 2^TU cWT H {e W B T_ ir 

(2) 

T{uj) is the transmission function, which in the 
ideal ballistic limit equals the number of phonon 
subbands, Nb(uj), at frequency ui (see Fig. [5J. In 
the low temperature limit we can approximate 
T(lo) = T(Q) = 4, since there are four acous- 
tic modes at hw — > |39j . The remaining in- 
tegral in Eq.© can be calculated analytically to 
yield the universal thermal conductance quantum 
G Q (T)^A(^P B T/3h). 

Figure [1] shows the thermal conductance, G, vs. 
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temperature for a 1.2 nm diameter (110) wire and 
a 1.0 nm diameter (100) wire. The results ob- 
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Figure 1. Thermal conductance, vs. tempera- 
ture for a 1.2 nm diameter (110) wire and a 1.0 nm 
diameter (100) wire. Solid lines: TEP without H- 
passivation. Dotted lines: TEP with H-passivation. 
Dashed lined: DFT. 

tained with the TEP model agree quite well with 
the microscopic DFT calculations. At tempera- 
tures T < 200 K, the relative error is less than 5%, 
and at T = 500 K the difference is within 10%. 
Also, observe that including hydrogen in the em- 
pirical potential calculations leads to a maximum 
deviation of 3% (at T = 50 K, (100) wire) as com- 
pared to the pure Si calculations. Leaving out the 
hydrogen in the TEP calculations is therefore jus- 
tified. We take these results as a validation of the 
relatively simple and fast empirical potential calcu- 
lations on pure Si wires, and proceed to study much 
larger wires - a study that would be very time con- 
suming using DFT calculations. In the rest of this 
paper, we study the diameter dependence of (100), 
(110) and (111) oriented SiNWs (see FigureEj), us- 
ing the TEP approach without H-passivation. We 
note that the same approximation was recently ap- 
plied in Refs. [HHD]. 

Figure [2] shows the atomic structure in the pe- 
riodically repeated unit cell (along the z-axis) of 
2 nm diameter wires with the three orientations. 
The top row shows a cross section of the wires while 
in the middle row we have shown side views of the 
unit cells. In the bottom panel the full phonon 
band structures for the three wires are shown. For 
each wire there are 3Na subbands, with Na being 
the number of atoms in the unit cell. Note that 
the different wire orientations have different unit 
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Figure 2. Cross sections (top) and side views 
(middle) together with the full band structures 
(bottom) of ~2 nm diameter wires with orienta- 
tions (100) (left), (110) (middle), and (111) (right). 
The wire is always oriented along the z— axis. The 
phonon wave vectors, q are all in the respective wire 
directions and are show in units of the reciprocal 
unit cell lengths, with ainrj = 5.4 A, ano = 3.8 A, 
and am = 9.4 A. 



cell lengths, ahki, and thus also different number 
of atoms in the unit cell. At low energies there 
are only four acoustic modes: the dilatational and 
the torsional modes with linear dispersions and two 
flexural modes with quadratic dispersions [39] . The 
many, relatively flat subbands starting right above 
the acoustic modes are denoted shear modes and 
originate from the bulk acoustic bands that are 
folded into the wire axis. The many closely ly- 
ing bands around E p h = 63 meV correspond to the 
bulk LO phonon energy. 

An important qualitative feature is that the 
bands in the (110) wire generally have a larger 
slope (velocity) than the the other wire orienta- 
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tions, where there are many flat bands. The same 
trend is also seen for other wire diameters. Due 
to the larger slope there are more bands at a 
given energy for the (110) wires, and one expects 



T(no) > T { 
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conductance for the (110) wire. 



leading to a larger thermal 
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Figure 3. Transmission function for wires 
with different orientations. The area of the wires 
are Aioo = 8.9 nm 2 , Ano = 11.2 nm 2 , and 
= 9.0 nm 2 . The three dashed curves show 
the function f(uj,T) (see text) at temperatures 
T = 5, 20, 100 K (from left to right). 

This anisotropy is further illustrated in Fig. [3] 
showing the transmission function for three differ- 
ent wires orientations, all having diameter ~ 4nm. 
While the (100) and (111) wires have approxi- 
mately the same transmission, the (110) wire has 
a significantly larger transmission at most energies 
above 10 meV. The three black dashed curves show 



the function f(to,T) = uj 
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appearing in eqOat temperatures T = 5, 20, 100 K. 
We observe that at T < 5 K, all wires have T(lo) = 
4, and display the quantized thermal conductance, 
G(T <5K)=G Q (T). 

The curve f(uj,T = 20 K) suggests that the ther- 
mal conductance of the three wires should be ap- 
proximately the same up to temperatures T k, 
20 K. At higher temperatures, the (110) wires will 
have a higher conductance. 

Figure |4] shows the thermal conductance per unit 
area, G/A, vs. cross sectional area, A, at differ- 
ent temperatures 0D]. The different symbols cor- 
respond to different wire orientations and different 
cross sectional shapes of the wires. At low tem- 
peratures (5 K, left) all wires have the same 1/A 
dependence and coincide with the analytical result 
Gq(T)/A shown as the dashed line in the graphs. 
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Figure 4. Thermal conductance per unit area, 
G/A, vs. cross sectional area at temperature T — 
5K (left), T = 20 K (middle), T = 300 K (right). 
The different symbols correspond to different wire 
orientations as shown in the legend. The dashed 
lines show Gq(T)/A. Note the different scale in 
the three plots. 



At T = 20 K (middle) , the area dependence of the 
three wires is still approximately the same as ex- 
pected from Fig. [3] Note that the thinnest wires 
still follow the quantum line, Gq(T)/A. This is 
because confinement lifts the bottom of the shear 
modes (see Fig. [5J and increases the energy range, 
where only the four acoustic modes exist. 

At T = 300 K (right) G/A is almost area inde- 
pendent. The room temperature area normalized 
conductances of the (110) oriented wires are ap- 
proximately 50% and 75% larger than the (100) and 
(111) wires, respectively. The small fluctuations in 
G/A for a given orientation are due to variations in 
the surface structure of the wires at different diam- 
eters. Notice also that for the (100) orientation, the 
square shaped wires have a larger conductance than 
the circular shaped wires. Again, we attribute this 
to the details in the surface structure. The room 
temperature conductance per unit area of the (100) 
and (111) wires are in fairly good agreement with 
values reported in Ref. [T2] and [13], respectively. 

Experimental [3] and theoretical [TU] works have 
shown that the thermal conductivity decreases with 
decreasing diameter, seemingly contradicting the 
results shown in Fig. 01 However, we have not in- 
cluded any scattering (surface roughnes, impurity, 
phonon-phonon, phonon-electron) in our calcula- 
tions, and our results thus represents an ideal up- 
per limit for the thermal conductance. Indeed, if 
we include scattering by vacancies we observe a de- 
creasing conductance per unit area with decreasing 
diameter |36j . Quantitative comparison with ex- 
periments should only be carried out in the fully 
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ballistic regime for short wires at relatively low 
temperatures where anharmonic effects are unim- 
portant, and for wires with very smooth surfaces. 

The anisotropies can also be obtained from the 
bulk band structure. We construct an approximate 
wire band structure by projecting the bulk subband 
surfaces onto lines, qn, along the (100), (110), and 
(111) directions. Figure \E\ (top) shows the orien- 
tation of the basic tetrahedrals together with the 
direction of the q\\ lines. The lines are equidistantly 
spaced on a fine square grid in the plane perpendic- 
ular to the wire direction, and are restricted to the 
first Brillouin zone. A total number of N lines gives 
67V projected subbands, since there are two atoms 
in the bulk primitive unit cell, and three bands per 
atom. From the projected band structure we obtain 
the transmission by counting the number of bands 
at a given energy. Each of the lines ranges from 
<7H = [0, ir/ahki], with dhkl being the unit cell length 
of a wire in the (hkl) direction. For a given cross 
sectional area, A, the different wire orientations will 
have different numbers of atoms in the unit cell and 
thus also different number of bands (see Figure [5]). 
To account for that, we scale the transmissions by 
a factor 3Na/6N, where Na = pAahki is the num- 
ber of atoms in a wire unit cell, and p = 0.05 A -3 
is the atomic density. 

Figure O (bottom panel) shows the transmis- 
sion functions obtained from the bulk band struc- 
ture. The large peak in transmission around 
hoj = 20meV for the (110) wire seen in Fig. [3] 
is reproduced from the bulk bands with the ac- 
tual transmission values being in close agree- 
ment. The room temperature thermal conduc- 
tances per unit area, G/A, obtained from Fig.[5]are 
G (100) /A = 0.59, G (ll0) /A = 0.77, and G {nl) /A = 
0.49 nW/(K ran 2 ), in qualitative agreement with 
the values in Fig. [4] Note, however, at low tem- 
perature the results are wrong, since the true wires 
have four acoustic modes, while the projected bulk 
band structures only have three. At higher temper- 
atures, this error is less important. 

Evidently, it is primarily the two TA subbands 
that are responsible for the large anisotropy. This 
is not surprising since the bulk band structure in 
the (110) direction has two non-degenerate TA sub- 
bands while in the (100) and (111) directions thcR-eferences and Notes 
TA subbands are degenerate. Thus, importantly, 

the anisotropy is not merely a property of the spe- (1) Hochbaum, A.; Chen, R.: 
cific silicon nanowires studied here, but can also Liang, W.; Garnett, E.; 



Figure 5. Top panel: Orientation of the ba- 
sic tetrahedrals for the (100), (110), and (111) di- 
rections. The arrow indicate the wire direction as 
well as the q\\ lines. Bottom panel: Transmission 
functions obtained from the projected bulk band 
structure. The different curves show the contribu- 
tion from the two transverse acoustic (TA) modes 
(dotted red), the single longitudinal acoustic (LA) 
mode (solid blue) and the three optical (LO+TO) 
modes (dotted black). The results corresponds 
to square shaped wires with cross sectional area, 
A = 9nm 2 . 



be expected to take place in other materials with 
the diamond structure. Indeed, projection of the 
germanium bulk band structure as described above 
yields similar results. 

In conclusion, we have shown that the ballis- 
tic thermal conductance in silicon nanowires is 
strongly anisotropic with wires oriented along the 
(110) direction having 50% and 75% larger conduc- 
tance per unit area than wires oriented along the 
(100) and (111) directions, respectively. We believe 
that the presented results will be relevant for future 
phonon-cngineering of nanowire devices. 
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